Magnetic Properties of the Superconducting State of Iron-Based Superconductors 
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We show that features of the dynamical spin susceptibility can unambiguously distinguish between 
different pairing symmetries of the superconducting state in iron-based superconductors. A magnetic 
resonance is a unique feature for the extend s^2y2-wave cos kx cos ky pairing symmetry. It is present 
in the pure superconducting (SC) state, but weakens in the mixed SC and magnetically ordered 
state. We calculate the the RPA correction to the NMR spin relaxation rate l/Ti and the Knight 
shift in the above states and show a good agreement with experimental results. Moreover, we argue 
that the energy dispersion of the magnetic resonance along c axis observed in neutron scattering 
experiments is also an indirect evidence supporting the s,^2y2 pairing symmetry. 
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Introduction - The recently discovered iron-based 
high-temperature superconductors[ll,[3,[l,|3,[3| have very 
intriguing magnetic properties @ which may hold the key 
to understanding their superconducting pairing mecha- 
nism. Theoretically, many possible gap pairing symme- 
tries have been proposed. Due to the proximity of the 
superconducting state to a coUinear antiferromagnetic 
state, as well as due to relatively weak phonon interac- 
tions, a magnetism-based superconducting mechanism is 
favored. Falsifying or verifying the prediction of pairing 
symmetry will be a significant step in understanding the 
nature of the magnetically driven SC pairing mechanism. 

Both weak and strong coupling approaches suggest 
that an extended s-wave pairing symmetry is favored 
0, H, Q. Based on a weak-coupling approach, an s- 
wave, (so called s^) stateQ in which the sign of order 
parameters changes between hole and electron pockets, 
is argued to be favored for repulsive inter-band interac- 
tions. However, the weak-coupling approach does not 
specify the exact form of order parameters, and is de- 
pendent on the degree of nesting between the electron 
and hole pockets. In the strong coupling approach, in a 
recent paper we showed that the pairing symmetry 
is determined mainly by the next nearest neighbor anti- 
ferromagnetic exchange coupling J2 [ill Il2 | and has 
an explicit form in momentum space, cos{kx) ■ cos(fcy). 
This result is model-independent, as long as the dom- 
inating interaction is next-nearest neighbor J2 and the 
Fermi surfaces are located close to the F and M points 
in the Brillouin zone. The cos{kx) ■ cos{ky) changes sign 
between the electron and hole pockets in the Brillouin 
zone. In this sense, it resembles the order parameter, s^, 
proposed through weak-coupling general arguments [1] . 

The magnitudes of superconducting gaps measured by 
ARPES at different Fermi surfaces are in a good agree- 
ment with the simple \cos{kx)cos{ky)\ form 13, 3 ist - 
but directly probing the sign change from electron to 
hole pockets still remains a big experimental challenge, 
despite several theoretical proposals [H, [l3, [ll [13] . For- 
tunately, the magnetic properties in the superconduct- 
ing state may provide us with indirect evidence of the 



pairing symmetry. Maier and Scalapino[20j first pointed 
out that, in the extended s-wave SC states, there is a 
strong coherent peak in the dynamic spin susceptibility 
x(Q,a;) at an energy hiu below the two-gap value and 
at a wavevector Q — (0, tt) identical to the magnetic or- 
dering wavevector in parent compounds. Recent neutron 
scattering experiments [111 [i^ [ij] have observed a mag- 
netic resonance in the superconducting states similar to 
the magnetic resonance observed in cuprates. This re- 
sult provides an indirect evidence of the extended s-wave 
pairing if the coherent peak is a unique feature of this 
pairing symmetry. Neutron scattering experiments [2^ 
have also observed a significant dispersion of resonance 
peak along c axis. 

Complementing neutron scattering, NMR is an impor- 
tant probe of magnetic properties. The experimental 
results of the NMR spin relaxation rate 1/Ti in these 
superconductors have been a challenge for theories pre- 
dicting the extended s-wave pairing symmetry due to 
the absence of coherent peak and the near cubic power 
law dependence on temperature 24, 25, [2^, 27 1. Sev- 
eral works [2I, [2^, [s^l have addressed this problem and 
pointed out that for the cos(fc2:) •cos(fc3^) order parameter, 
the inter-band contribution to the NMR spin relaxation 
rate does not exhibit a coherence peak. As such, the 
experimental result is not inconsistent with the pairing 
symmetry if the inter-band contribution is much larger 
than the intra-band contribution or if the samples are 
strongly disordered. However, in a simple mean-field 
state, the intra-band contribution is always larger or com- 
parable to the inter-band contribution and the spin re- 
laxation rate still exhibits an enhancement^right below 
Tc owing to its fully gapped s-wave nature j28j. 

In this paper, we perform a detailed calculation of the 
magnetic response in the superconducting state of iron- 
based superconductors. Within a two orbital model and 
within the RPA approximation we address several exper- 
imentally related issues: (1) We demonstrate that the 
coherent resonance peak is an unique feature of the ex- 
tend s-wave pairing as suggested in ref. 2^, 3lj . We show 
that the different pairing symmetries in iron-based super- 
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conductors are distinguishable by the distinct features of 
their corresponding dynamical spin susceptibilities; (2) 
We investigate the magnetic susceptibility in the mixed 
SC and SDW state and find that the coherent peak is 
strongly weakened in the mixed state. Many experiments 
suggest a possible coexistence of SC and magnetic order- 
ing in underdoped materials [H, [H, 34 1, and the cal- 



culated doping dependence of the resonance peak could 
be checked in future experiments; (3) We show that by 
considering the RPA correction to the mean-field BCS 
state, both the 1/Ti and the Knight shift show a good 
agreement with experimental results even when the intra- 
band contributions are included. The RPA correction of 
the mean-field BCS state can dramatically enhance the 
inter-band contribution to the NMR spin relaxation time 
1 /Ti ; (4) We show that the dispersion of the magnetic 
resonance along the c axis observed in neutron scatter- 
ing experiments [23} can also be interpreted as indirect 
evidence of the Sx2y2 pairing symmetry. The magnetic 
exchange coupling along c-axis has two effects on the 
resonance: (1) modifying the SC gap along c-axis and 
(2) causing new RPA corrections. Although the second 
effect produces the right trend of the dispersion, it is too 
small to explain the dispersion observed in experiments. 
Therefore, it is the first effect that causes the disper- 
sion. Moreover, in order to explain the experimentally 
observed dispersion of the resonance peak towards lower 
energy from Qz = to Qz = tt, the superconducting or- 
der parameters at electron and hole pockets must have 
opposite values. 

Model- The Fe ions form a square lattice in the FeAs 
layer of LaFeAsO system with As ions sitting in the 
center of each square plaqucttc of the Fc lattice and 
displaced above and below the Fc plane. The crystal- 
lographic unit cell includes two Fe and two As ions. 
We adopt the two orbitals per site model proposed in 
Ref . [sit . capturing the degeneracy of the d^z and dyz 
orbitals on the Fe atoms and the general location of 
the Fermi surfaces. The resulting Fermi surface con- 
sists of two hole pockets, centered around the F = (0,0) 
and (±7r, ±7r), and two electron pockets, centered around 
(±7r, 0), (0, ±7r) in the unfolded Brillouin Zone (the hole 
pocket at (±7r, ±7r) is an artifact of the two-band model, 
all the LDA calculations as well as more complicated 5- 
band modelsfsi^l have this pocket at the F point). The 
single electron Hamiltonian is written as 
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(1) 



where -0^^ = {c\^^,c\^^) is the creation operator with 
spin-(T in the orbitals (1, 2) = {d^z^ dyz), fJ. is the chemical 



potential. The electronic dispersions are given by 

ex(k) = — 2ti cos kx ~ 2t2 cos ky — At^ cos kx cos ky, 
(k) = — 2<i cos ky — 2t2 cos k^ — 4t3 cos k^ cos ky, 
£2,1, (k) = — 4t4sinA:a; sinfcj,. 

(2) 

Choosing ti — —1, t2 — 1.3, and ~ t4 — —0.85 in 
electron volts, the eigenvalues of Eq. (P), 



(3) 



yield the two Fermi pockets aforementioned. For the half- 
filed system, corresponding to two electrons per site, fi — 
1.54. 

In the following we shall consider the singlet pair- 
ing between electrons within each orbital based on t — 
Ji — J2 model proposed in ref.Q- The symmetry of the 
superconducting order parameter A(k) has two possi- 
ble d-wave types, d^^-y^ = Ao(cosfc2: — cosky)/2 and 
dxy = Aq sin fca; sin fcj, , and two s-wave types, s^^+y^ = 
Ao(cosfcx + cosky)/2 and Sx2y2 — Aq cos fc^ cos fcj,. We 
neglect dxy and the inter-orbital pairings: ref. 0] showed 
they are negligible for the case oit — Ji — J2 model. The 
effective mean-field Hamiltonian is then given by 



i?MF = ^«'I/l(k)«'k, 



(4) 
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component spinor = (cj ,^.1., ci _k,i, 4 ^ C2,-ka)- 
Due to the C4 symmetry of the underlying lattice, 
Ai(k) ~ A2(k), except for the dx^-y-i case, where 
Ai(k) — — A2(k). The dx^^y^ and Sx^j^yi pairing sym- 
metries are nodal and Sx^yi is nodeless for any small dop- 
ing parameter. Note that 52,2^2 exhibits a sign change 
between the hole and the electron pockets, while Sx^j^yi 
does not. 

First we consider the one-loop contribution to the spin 
susceptibility that include the intra-band and inter-band 
contributions as 

xt^^M = dr{TSr{q,r)Sf{-q,0))e'-- (6) 

where Sf{q,T) = Efe cj,_k+q,^ Wc'Mi.Ca,k,i.(T), and a, (3 
refer to different orbital indices. The total susceptibility 
is given by Xo (q>«'^) = T,ai3Xif{(l,'i^)- Since x"^" = 



for the singlet pairing, we calculate 



(X"" + X''^)/2 = X^ 
the susceptibility, 

Xo" (q, [^('^ + ^' + »^)G(k, lu;^)] , 

k,n 

(7) 
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where G{]s.,iuj) = {iujl — Hmf)~^- Introducing the band 
operators, $k, such that = i7(k)$k, the Haniiltonian, 
HyiFj is diagonahzed with the eigenvalues, i?i(k), and 
the eigenvectors, Af^''^(k, q), and the spin susceptibihty 
becomes 



Xo'(q,«w) = 



-i^M^(k,q)x 

k,a,/3 

/(£:,(k))-/(£:,(k + q)) 



iLU — £j (k) + £i(k + q) 



(8) 



where f{Si{k)) = 1/(1 + exp[fi(k)/T]). For s-wave sym- 
metries, Ai ~ A2, after the analytical continuation we 
obtain: 



xo^(q,'^) 



E 



A(k)A(k + q) + g^(k)^„(k + q) 



r 



/(g,.(k))-/(-g,(k + q)) 
LO + ir] - £f_i{k) - £^(k + q) ' 



(9) 



where £±{k) ^ ^ El{k) + A^{k). 

Within RPA, the spin susceptibihty XRPA(q, is 
given in matrix form by: 

XRPA(q, itj) = [I - rxo(q, icj)]" Vo(q, i^j), (10) 

where I is a unit matrix and the vertex T — Ul with 
the value of U chosen in the paramagnetic phase, xo 
above is written in 2 x 2 matrix form whose entries con- 
tain the orbital contributions of Eq. ([7]). Note that the 
RPA enhancement of the spin fluctuations is determined 
by the det|I — rxo(q, In the superconducting state 
the quasiparticles at the hole and electron Fermi surfaces 
are connected by the Q = (7r,0), (0,7r) wave vector. In 
the unfolded BZ, the 3^,2^2 order parameter satisfies the 
condition A(k) = — A(k-I-Q). The imaginary part of 
the inter-band bare susceptibility is zero for small fre- 
quencies due to the opening of the gap and experiences a 
discontinuous jump at uJc = min(|A(k)| -f |A(k -I- Q)|, k). 

Comparison of spin susceptibility in different pairing 
symmetry states: Fig. [1] presents the imaginary part of 
the RPA spin susceptibility ImxRPA(Q, "j^) for various 
superconducting order parameters on (a) the electron- 
doped side and (b) the hole-doped side. We observe 
that only the s^^y^ pairing symmetry gives peaks un- 
suppressed by RPA. At both electron and hole dop- 
ings, there are two pronounced peaks for this pair- 
ing symmetry. While the number of peaks is model- 
dependent and in this case characteristic to the two- 
band model used, the high intensity of the peaks is a 
generic feature stemming from the coherence factor (1 — 
A(k)A(k+Q)/£:+(k)£:_(k+Q))/2, where Q = (7r,0). 
This coherence factor clearly becomes large if A(k) = 
— A(k+Q), a condition only met by the 52,2^^2 pairing. 
The coherence factor is weighted by the spectra weight 
S{iLj — ^^+(k) — £_(k+Q)) and integrated over k to ob- 
tain the final Im{x)- For the Sx2y2 symmetry pairing, 
the energy profile (after the SC gap opening) is almost 
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FIG. 1: The imaginary part of the RPA spin susceptibil- 
ity xrpa at Q — {n,0) for various superconducting pairing 
symmetries, (a) hole-doped side {fi — 1.4) and (b) electron- 
doped side (/i = 1.6). The blue line represents the result 
for s^2j^2 ~ cos kx cos ky, the green dashed line represents the 
Sj.2^y2 ~ cosfca; + cos ky, and the red dashed line represents 
the dx2_y2 symmetry. We used Aq — 0.3 and U = 0.6ti. 



isotropic and flat around the F and X points. This means 
that at Lu UJc, almost all the k's around the electron and 
hole FSs can contribute to the coherent peak, making its 
intensity significantly higher than the other pairing sym- 
metries. For d^2_y2 {s^2_^_y2) pairing symmetries, there 
are four nodes on around F {X) point, making the en- 
ergy profile around the F (X) point considerably angle- 
dependent. Therefore it is rather difficult to find a con- 
stant UJc such that most k's on the FSs can contribute to 
the sum. Apart from the number of peaks, the previous 
discussion is generic and does not depend on the particu- 
lar model used. Within the two-band model, the number 
of peaks for the s,j.2y2 pairing symmetry is two. This is 
due to the FS topology given by the two band model. 
As shown in FIGl2]the two resonance peaks corresponds 
to the nesting of two pairs of electron-hole pockets: the 
nesting between F hole pocket and X electron pocket 
and the nesting between M hole pocket and Y electron 
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FIG. 2: The schematic of the FSs not far away from half- 
fiUing. The hole and electron pockets are connected by the 
vector Q. 



pocket. When doping is small, we can roughly express 
the two resonance frequencies as 



Wros,2 



|A(k^) 
|A(k2^) 



|A(M)|, 
|A(k^)I. 



According to the pairing gap formula cos(fca;) cos(fcj,), the 
larger FS gives a smaller gap magnitude. Therefore we 
have Wros.i > '.i^rcs,2, bccausc the FS around F is larger 
than that around M. For the d^^^y^ pairing, the two 
hole pockets have nodes and the minimum gaps A(k'')'s 
on the hole surface are hence much smaller than the elec- 
tron surface gaps A(k^)'s; for the s^^+i/^ pairing, the two 
electron pockets have nodes and the gap values on the 
electron pockets are much smaller. Hence, the peaks for 
these two pairing symmetries shift to the lower energy 
than to those of the Sr^2y2 symmetry. 

The band structure of the current two band model 
as the function of doping will not exactly represent the 
real materials. However, the differences of the resonance 
peaks between electron-doped and hole-doped strongly 
suggests that the resonance peak is tied to the combi- 
nation of the opposite signs of SC order parameters and 
the close matching of Fermi surfaces between electron 
and hole pockets. The resonance is stronger for bet- 
ter nested surfaces connected by Q vectors. This pro- 
vides a direct explanation of the orbital selective electron- 
mode coupling observed in Angle Resolved Photoemis- 



sion Experiments 37 1 . 

Magnetic resonance in the mixed SC and SD W states 
with Sx2y2 pairing symmetry Recently, several experi- 
ments indicate the coexistence of SC and SDW ordering 
in underdoped materials 3^ 33, 34|. Therefore, it is of 
great interest to investigate whether the magnetic reso- 
nance survives in the mixed state. Treating the SDW 
order in mean-field, we now have off-diagonal terms be- 
tween k and k -|- Qsdw, where Qsdw = (tt, 0). In the 
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FIG. 3: The imaginary part of the RPA susceptibility in the 
(a)hole-doped and (b)electron-doped side. We used the same 
Ao and the chemical potentials with Fig. [l] 



(11) mean-field approximation, the spin exchange interaction. 



ffcx= ^(q)S(q) • S(-q), (12) 

can be decoupled as 

= E ^(q)(s(q)) • s(-q) + h.c, (13) 

q 

where J(q) = Ji(cos(gj;) -I- cos((jy)) + 2 J2 cos(g^) cos(gy). 
Substituting the realistic magnetic structure, (S(q)) — 
(0, 0, Af)(5q.Q3„„, , the mean-field interaction reduces to 

^^eT = W^E4k,.<^'c,.k+Q,,„..' +h.C., (14) 

where W = J{QIsdw)M is a control parameter chosen 
between ~ Aq in our calculation. Therefore we can 
write the full mean-field Hamiltonian as 



(15) 



where 



J 



/i(k) = /i(k) (/ -h T^)/2 + /i(k q,sDw) r^)/2 + WI®I<^T^, 



(16) 



(17) 
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with / a 2x2 identity matrix and /i(k) a 4x4 matrix 
given in Eq.®. 

Fig. [3] presents the RPA total spin susceptibihty in the 
mixed state of SDW coexisting with s^2y2 order in the (a) 
hole-doped and (b) electron-doped side. The pronounced 
resonance peaks observed in the pure Sx2y2 superconduct- 
ing state are reduced by the presence of the SDW. The 
SDW order doubles the unit cell and reduces the orig- 
inal Brillouin zone which causes the weight of the spin 
susceptibility at (0,7r) to spread out. 

NMR spin-relaxation rate and Knight shift As dis- 
cussed earlier in the manuscript, the experimental ab- 
sence of a coherent peak and the presence of near cubic 
power law dependence on temperature of the NMR s pin 
relaxation rate 1 /Ti in these superconductors 24 , 2^ 2^, 
271 have been a challenge for theories predicting the ex- 
tended s-wave pairing symmetry. In a previous paper (28j|. 
two of us also calculated the NMR spin relaxation rate 
l/Ti of the bare superconductor and found that it fac- 
torizes into inter- and intra-band contributions. While, 
for the cos(fca;) • cos{ky) order parameter, the inter-band 
contribution to the NMR spin relaxation rate does not 
exhibit a coherence peak, the intra-band contribution is 
larger than the inter-band contribution and still exhibits 
an enhancement right below Tc owing to its fully gapped 
s-wave nature. Adding the two contributions we find 
that, although the coherence peak for cos{kx) ■ cos{ky) 
is smaller than that for a sign-preserving gap (such as, 
for example, | cos(fcx) • cos(fcy)|), it is still present due to 
the intra-band contribution. The coherence peak can be 
strongly reduced if the intra-band scattering is stronger 
than inter-band scattering. Here we also show that the 
NMR experimental results are consistent with the ex- 
tended s-wave pairing symmetry if the simple RPA cor- 
rection over the mean-field bare SC state is considered as 
argued by 29]. 

The spin-relaxation rate at a temperature T is defined 



TiT 27r i^o^o cuq 



- - - without RPA 
within RPA 




FIG. 4: (a)The NMR spin-relaxation time l/Ti as a function 
of temperature T in unit of Tc for various values of F, which is 
considered as the phenomenological disorder parameter. We 
used Ao/2 = Tc with Aq = 0.3. (b)The effect of the spin 
fluctuation within RPA susceptibility. We choose the on-site 
Coulomb repulsion U — 0.5ti such that the system remains 
in the paramagnetic state. Inset: the log-log plot of the RPA 
susceptibility with line of T^ showing the temperature depen- 
dence of xrpa just below Tc. 



normaf 

cos kx + cos ky 




0.2 0.4 0.6 0.8 1 

T/Te 

FIG. 5: Temperature dependence of the Knight shift. The 
solid green line is a plot with a sign-reversed extended s- 
wave, cos kx cos ky, and the dashed red line is one with sign- 
preserved extended s-wave, cosfe^ + cos ky. Inset shows the 
difference between two order parameters just below Tc. 



where ^(q) is a structure factor: For F it is roughly 
isotropic, while for As, A{q) = cosqx/2coaqy/2 (neglect- 
ing the fact that the position of As is out of plane of Fe's). 
In the following we take the structure factor uniform, i.e., 
A(q) = 1. 

In Fig. m we present the NMR spin-relaxation time as 
a function of the temperature(in unit of T^). The coher- 
ence peak is present in the bare Sx2y2 pairing SC state 
if the impurity effect is weak, as found in [28[. The co- 
herence peak is reduced by the effects of disorder and 
the RPA correction. The disorder is taken in account by 
broadening the imaginary part of frequency, F. Fig.[4^a) 
shows that the coherence peak at Tc is reduced as dis- 
order increases, F — Tc/10 ^ Tc/3. Fig. [DJb) presents 



the effect of the RPA correction. A small interaction 
U ^ 0.5ti is enough to suppress the coherence peak. 

We also calculate the Knight shift. Kg = x(ct ^ 
0,uj = 0). Fig. O plots the Knight shift value as a 
function of temperature in the different s-wave pair- 
ing symmetry states: a sign-reversed extended s-wave 
{Sx2y2 ~ COS k^ cos ky) and an extended s-wave (sj,2_|_j,2 ~ 
cosfcj; -I- COS ky). The main difference between these two 
ordering symmetry states is that the former one is fully 
gapped while the second one has nodes on the electron 
pockets. At low temperatures, KgiT) has a power-law 
temperature dependence for the second case. The fact 
that Ks{T) is a concave function around T = 0.4Tc in- 
dicates the presence of two gap values, as observed in 
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FIG. 6: (a) The RPA spin susceptibilities in the electron- 
doped s^2y2 superconducting state(/i = 1.65|ti|) with a var- 
ious ?7int(q), and (b) the dispersion of a function of 
t/i„t(q). 



experiments [2J| . 

The effect of magnetic exchange coupling along c-axis: 
Now we turn to the effects of the spin fluctuations be- 
tween the FeAs layers. For the 122 series of iron-based 
superconductors, the coupling between FeAs layers along 
c-axis is very important. In the parent compounds, 
the magnetic exchange coupHng SJz along c-axis de- 
termined by neutron scattering experiments is around 6 
meV - around one fifth of the in-plane magnetic exchange 
coupHngjsll; this suggests a strong three-dimensional 
electronic structure. The measurement of spin excitation 
resonance peaks in BaFe1.9Nio.1As2 shows they are dis- 
persive along c-axis direction, and the resonance peak at 
3D AF ordering wavevector (tt, 0, it) is hoj — 7meV below 
Tc different from one at Qi = (tt, 0, 0), fiw = 9.1meVl23i]. 

The Hamiltonian for the antiferromagnetic exchange 
coupling in c-axis can be generally written as 

Hint = J.^S"(r).S"(r±z) 

+ JiY. S"(r).S'^(r±z), (19) 



where Jz is a spin coupling constant for the intra-orbital 
and J'z for inter-orbital between the adjacent FeAs layers. 

In general, there are two possible effects generated by 
this new spin coupling. First, it could produce the varia- 
tion of the superconducting gap as a function of momen- 
tum wave- vector along c-axis. This effect has been used 
to explain the dispersion of magnetic resonance [3] . Here 
we point out that the explanation is only consistent with 
the extended Sx2y2 wave. Let's revisit the analysis given 
in ref-lil]. Let Ag and Ag denote the superconducting 
gaps on the hole and electron pockets respectively (in a 
pure two dimensional model). By considering the antifer- 
romagnetic coupling between layers, the gap functions, at 
the mean- field level, are naturally modified to Ae{kz) = 
AO + Scos{kz) and Ah{kz) = ^1 + 5cos{kz). For a S"^ 
pairing symmetry, Ag ~ — Aq ^ — Aq. Therefore the 
dispersion of the resonance along c-axis is roughly deter- 
mined by 



hu{qz) - mm{{\Aeikz)\ + \Ahikz + qz)\),kz 
~ 2Ao-2<5|sin(|)| 



(20) 



In the above analysis, if we assume the sign of gap does 
not change from electron to hole pockets, it is easy to 
see that the dispersion of the resonance in Eg 1^ changes 
to huj{qz) ~ 2Ao -I- 2(5|sin(^)|, which will result in an 
opposite dispersion of the resonance energy than that 
reported in the experiments: larger at the wavevector 
(tt, 0,7r) than at the wavevector (7r,0, 0), obviously con- 
tradicting the experimental results. 

There is also a second effect due to the magnetic ex- 
change coupling along c-axis, namely, a simple RPA cor- 
rection to spin susceptibility due to the presence of the 
exchange coupling. It is easy to show that in the presence 
of above exchange coupling, the RPA spin susceptibility 
is modified to 



X^P^(q,c.) = ^[xo(q,c.)(l - V{cDxo(.ci,Lu))-%0, 

a/3 



with the vertex 



V(n) = ( -J'zCOSqz 

"^^"^^ \ -J', COS qz C/int(q) 



(21) 
(22) 



where C/int (q) = U — Jz cos qz with U a onsite Coulomb 
repulsion. We calculate xrpa(w, q) in the electron-doped 
superconducting state(^ = 1.65|ti|) with Sx2y2 pairing 
symmetry, A(k) — Aq cos fc^, cos fcy. J^, the inter-orbital 
coupling has very little effect on the positions of Wros- 
Therefore, we only need to consider the intra-orbital spin 
coupling between the layers, Jz- The effect of Jz simply 
creates an effective Uint{q) which is larger at g = (tt, 0, tt) 
than at g = (vr, 0, tt). Fig. Mj^) displays the dispersion of 
ujj-cs as a function of C/int(q)- It is clear that the resonance 
energy decreases as t/int increases. Therefore, in general. 
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the RPA correction would result in a dispersion of reso- 
nance peak that has a correct trend as the experimental 
results. However, it is also clear from FiglS]that Jj has 
to be comparable with U in order to create a visible dis- 
persion. Since Jz is really a small energy scale compared 
to both band width and interaction, the dispersion due 
to RPA correction from Jz is negligible. 

From above analysis, we can conclude that the dis- 
persion of the resonance peak along c-axis mainly stems 
from the modification the SC gap along c-axis. The fact 
that the dispersion of the resonance peak towards lower 
energy from Qz = to Qz = tt indicates that the super- 
conducting order parameters at electron and hole pockets 
must have opposite values. 

Conclusion We have investigated the magnetic prop- 
erties in the SC states of the iron-based superconductors 
based on a two-orbital model, identifying generic, model- 
independent, experimentally observable predictions. We 
have found the Sx2y2 pairing symmetry exhibits a strong 
magnetic resonance, which is absent from other pai ring 
symmetries (the similar result was also reached in |39j 
when this paper was prepared). We have predicted that 
the coexistence of SDW order with the SC order weak- 
ens the resonance peak. We have shown a good agree- 
ment with NMR experimental results, both the 1/Ti and 
the Knight shift, can be reached after considering the 
RPA correction from a bare BCS state with the extended 
Sx2y2 pairing symmetry. We have also explicitly studied 
the dispersion of magnetic resonance along c-axis [13 and 
have concluded that the dispersion is mainly caused by 
the modification of SC gap along c-axis and it indicates 
that the superconducting order parameters at electron 
and hole pockets must have opposite values. This ex- 
planation is also consistent with the three-dimensional 
nature of the 122 series of iron-pnicitides as shown in 
different experiments [H, |4^. Finally, we also want to 
point out that this explanation does not contradict the 
possible existence of gap nodes along c axis which may 
explain many nodal-like experimental results in the 122 
materials |4]|, since the resonance peak is determined 
by the summation of the absolute values of SC gaps at 
electron and hole pockets. 
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